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ABSTRACT 

We report the discovery of a radio counterpart to GRB 990123. In contrast 
to previous well-studied radio afterglows which rise to peak flux on a timescale 
of a week and then decay over several weeks to months, the radio emission from 
this GRB was clearly detected one day after the burst, after which it rapidly 
faded away. The simplest interpretation of this "radio flare" is that it arises 
from the reverse shock. In the framework of the afterglow models discussed to 
date, a forward shock origin for the flare is ruled out by our data. However, 
at late times, some radio afterglow emission (commensurate with the observed 
late-time optical emission, the optical afterglow) is expected from the forward 
shock. The relative faintness of the observed late-time radio emission provides 
an independent indication for a jet-like geometry in this GRB. We use the same 
radio observations to constrain two key parameters of the forward shock, peak 
flux and peak frequency, to within a factor of two. These values are inconsistent 
with the notion advocated by several authors that the prompt optical emission 
detected by ROTSE smoothly joins the optical afterglow emission. Finally, with 
hindsight we now recognize another such radio flare and this suggests that one 
out of eight GRBs has a detectable radio flare. This abundance coupled with 
the reverse shock interpretation suggests that the radio flare phenomenon has 
the potential to shed new light into the physics of reverse shocks in GRBs. 

Subject headings: gamma rays: bursts; radio continuum: general; shock waves; 
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1. Introduction 

BeppoSAX ushered in 1999 with the discovery of a super-bright 7-ray burst, 
GRB 990123. This GRB was intensively studied by many groups worldwide. The burst is 
notable for the richness of new results: the discovery of prompt optical emission by ROTSE 
(Akerlof et al. 1999), the discovery of the brightest optical afterglow to date and its redshift 
(Kulkarni et al. 1999), and a break in the optical afterglow light curve (Kulkarni et al. 
1999, Fruchter et al. 1999, Castro-Tirado et al. 1999). 

To date, as a result of our radio afterglow program at the Very Large Array (VLA)Q 
and the Australia Telescope Compact Array (ATCA) we have detected radio afterglow 
emission in nine GRBs (e.g. Frail et al. 1997b; Taylor et al. 1998). Observations of radio 
afterglow have resulted in the first direct proof for relativistic expansion in GRBs (Frail et 
al. 1997b) and robust estimates of the total burst energy (Frail, Waxman & Kulkarni 1999). 

The radio afterglow, when detected, exhibit the following pattern: the flux rises 
gradually to a maximum, persists for weeks to months before fading away. It is against this 
backdrop that we report, in this Letter, the discovery of a new phenomenon in the radio 
afterglow of GRBs - the discovery of a short lived (~30 hours) radio emission in GRB 
990123. This discovery, in a manner similar to the ROTSE discovery, adds yet another 
diagnostic - a probe of the reverse shock - to the study of GRBs. 



7 The VLA is operated by the National Radio Astronomy Observatory which is a facility 
of the National Science Foundation operated under a cooperative agreement by Associated 
Universities, Inc. 
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2. Observations 

OVRO Millimeter Array. Continuum observations in the 3-mm band were made with 
the Owens Valley Radio Observatory (OVRO) six element array on 1999 January 28.39 
UT. We obtained about 5 hours of on-source integration time in good 3-mm weather. Two 
1-GHz bandwidth continuum channels with central frequencies of 98.48 GHz and 101.48 
GHz were observed simultaneously in the upper and lower side bands. The quasar 3C345 
was used for phase calibration and 3C273 for absolute flux calibration; the observing 
procedure was otherwise similar to that described in Shepherd et al. (1998). 

OVRO 40-m Telescope. The position of the optical transient was observed with the 
newly commissioned 30-GHz receiver on the 40-m OVRO telescope. The receiver has two 
beams, each with a full- width at half maximum (FWHM) of 1.5 arcminute, separated by 
7.6 arcminute on the sky. The two beams were switched at 125 Hz and flux calibration 
was achieved by observations of 3C286. Pointing was monitored by frequent (20 minutes) 
observations of J1419+383 and J1613+342. 

JCMT Observations in the 1.35-mm and 850-/zm bands were made at the James Clerk 
Maxwell Telescope (JCMT)0 using the Submillimeter Common-User Bolometer Array 
(SCUBA; Holland et al. 1998). The atmospheric opacity in the 850 fim band was determined 
from skydip observations and that in the 1.35 mm band was inferred from skydips made at 
225 GHz using a tipping radiometer. The short term variability of the atmospheric emission 
was estimated from the outer ring of pixels and subtracted from the central (photometric) 
pixel. Frequent observations of a blazar established that the pointing of the detector was 

8 The JCMT is operated by The Joint Astronomy Centre on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, the Netherlands Organization for 
Scientific Research, and the National Research Council of Canada. 
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better than 2 arcsecond. Except for the observations on January 24, the rest of the data 
were taken under excellent sky conditions. 

VLA. All observations were obtained in two independent 50-MHz channels centered on 
8.43 GHz and 8.48 GHz. The interferometer phase was calibrated using the nearby point 
source J1545+478. The flux scale was tied to the source 3C286. Data calibration and 
imaging were carried out with the AIPS software package following standard practice. A 
log of the observations and a summary of the results can be found in Table p]. 



3. Results 

OVRO Millimeter Array. The source was not detected on January 28. The final 
image has an effective central frequency of 99.98 GHz, a synthesized beam of 1.25" x 0.83" 
(FWHM) at position angle 88.1°, and an rms noise level ~ 1.1 mJy/beam. At the position 
of the OT we find a peak brightness of — 1.44 mJy/beam. 

OVRO 40-m. Each observing run was about 2 hours long and the epochs of the four 
runs were (UT) January 24.5, February 5.5, February 6.5 and February 7.5. Bad weather 
prevented us from reaching thermal limits by almost a factor of 10. Averaging the four 
datasets we can place 3-er upper limits of 13 mJy (Band A, 26-28 GHz), 10 mJy (Band B, 
28-30 GHz), 12 mJy (Band C, 30-32 GHz) and 17 mJy (Band D, 32-34 GHz). 

JCMT. As can be seen from Table the source was not detected at any epoch save perhaps 
that of February 4.81 UT. On this epoch a source is seen with a signal-to-noise ratio 
(SNR) of 3.3, nominally a significant detection. However, single dish observations, unlike 
interferometric observations, are susceptible to several systematic errors. Our confidence in 
this detection would be strengthened if the if the source was well detected at least at one 
other epoch. 
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VLA. A strong source was detected in only one observation (January 24; see Table 0). 
This unresolved source, hereafter VLA J152530. 3+444559, is located at (epoch J2000) a = 
15 h 25 m 30.31* (+0.0P) 5 = +44 45'59.24" (±0.15). The coordinates of this source are in 
excellent agreement with the position of the optical transient (Kulkarni et al. 1999). As 
can be seen from Figure [l] the source rapidly fades and is only marginally detected (at the 
2-0" level) over the next few weeks. Other than the detection on January 24th, there was no 
reliable detection at any other frequency or at any other epoch including other radio efforts 
(Galama et al. 1999). To our knowledge, this phenomenon - a short-lived radio emission or 
flare - is a new phenomenon in radio afterglow studies. 

Given our claim of the discovery of a new phenomenon, it is important for us to 
demonstrate confidence in the existence of the source VLA J152530. 3+444559. To this 
end we make the following observations. (1) The source is seen in many different subsets 
of the data of January 24. (2) The source is seen in both (circular) hands of polarization 
(the fractional circular polarization is less than 37% at the 99.9% confidence level). (3) The 
source is separately detected in each of our 50-MHz wide channels. (4) Finally, within errors 
the source has the expected point-spread function. In view of these tests, we conclude that 
the source VLA J152530. 3+444559 is indeed a real source. 

Around the epoch of the radio discovery, both optical (Kulkarni et al. 1999; Castro- 
Tirado et al. 1999, Galama et al. 1999) and X-ray afterglow (Heise et al. 1999) were strong 
and fading. In order to compare the radio observations with those at other wavelengths we 
parameterize the temporal behavior of the 8.46-GHz flux by a two-component power-law 
model: flux at time t, f(t) = /*(t/£*) Qr ' d - Here, the time t is measured with reference to the 
GRB event and t* = 1.24 d is the epoch of the VLA detection of the radio emission; /* is 
the 8.46-GHz flux at that epoch and the subscript "r" ("d") stands for "rise" and refers to 
t < t* ("decay" with t > £*). Between January 23.63 and 24.65, we find the flux rose with 
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a r = 0.83 (mean estimator). From a Monte Carlo analysis we find 0.45 < a r < 1.85 (95% 
confidence level). For t >t* we find /* = 250 //Jy and = —1.3 (14 d.o.f., \ 2 = H-6). At 
the 95%-confidence level, a& < —0.9 and is < —0.8 at the 99.9%-confidence level. Clearly 
the source faded and we can, with great certainty, rule out light curves (for t > i*) with 
rising, constant or even slowly decaying emission. 

4. Origin of the Radio Flare 

The simplest possibility is that the radio flare was amplification due to interstellar 
scattering and scintillation (ISS). However, if this was a typical radio afterglow we 
expect the afterglow to grow stronger with time along with additional opportunities for 
scintillation. As can be seen from Figure [l] the source was never reliably detected after 
January 24 thereby rendering the ISS hypothesis unlikely. 

Sari Sz Piran (1999) suggest that both the prompt optical emission (Akerlof et al. 1999) 
and our radio flare arise from the reverse shock. In this model, the particles in the reverse 
shock cool adiabatically while there is a continual supply of newly energized particles in the 
forward shock. Consequently, with the passage of time, the emission from the reverse shock 
shifts rapidly to lower frequencies. At early times, self-absorption suppresses the radio flux 
and at late times adiabatic cooling results in weak emission. Due to these competing effects 
the emission from reverse shock is seen briefly - a radio flare. In contrast, Galama et al. 
(1999) have proposed that the radio flare arises in the forward shock i.e. the usual afterglow 
emission. We now investigate this proposal in some detail. 

At a given time, the broad-band afterglow spectrum can be described by three 
power laws and three characteristic frequencies (e.g. Sari, Piran & Narayan 1998): (z) 
/„ oc z/~( p_1 )/ 2 at high frequencies, v > v m \ (ii) f v oc v 2 at low (radio) frequencies, v < v a 
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with (Hi) f v oc z/ 1 / 3 bridging the two regimes. Here, p is the power law index of the shocked 
particles in the forward shock and is usually about 2.5. v m is the frequency at which the 
broad-band spectrum peaks (peak flux is indicated by f m ) and v a is the frequency at which 
the synchrotron optical depth is unity. 

As noted by Galama et al. (1999), the multi-wavelength data at the epoch of the VLA 
8.46-GHz detection (epoch, £*) can be fitted to an afterglow model with v m ~ 500 GHz and 
v a = 7.5 GHz. In this model, one expects the 8.46-GHz flux to initially rise oc t 1//2 and decay 
when u m falls below v a . This is expected to take place at epoch £*(500 GHz/V*) 2 / 3 ~ 19i*; 
here, = 8.46 GHz. This expectation is in clear contradiction with our measurements (§0). 

The rapid fall in the radio data thus requires that v m be comparable and preferably 
below v*. Galama et al. (1999) noting this requirement suggested a model with 
v m ~ v a ~ 30 GHz at epoch i*. In Figure ^| we display the spectrum with v a = v m = 39 
GHz at epoch £*; hereafter, we refer to this as model "Fq". The specific value of v m was 
chosen so the spectrum fits both the optical r-band flux and the VLA 8.46 GHz detection. 
Model Fg violates the 15 GHz and 86 GHz upper limits. Furthermore, in this model, we 
expect the radio flux to start decaying at epoch t*(^ m /^*) 2 / 3 or 3.5 d - in violation of the 
VLA observations. The true situation is even worse. When v m < u a , the latter will not 
remain time invariant but evolves as v a oc t _ ( 3p+2 ^ 2 ^ p+4 ' ) = t~°- 73 . This should cause the 
flux to increase faster as t 5//4 . Only when v a falls below v* will the 8.46-GHz emission start 
decaying and this will not happen until 8.2£* = 10 d. The above argument used asymptotic 
limits of the broad-band curve. In Figure [I] we display the evolution of the models with 
spectra shown in Figure using the usual prescription (e.g. Sari et al. 1998). The refined 
light curves, while smoother, do not agree with the observations. 

To conclude, in the framework of standard afterglow models (by which we mean models 
which assume a power law shocked particle spectrum and constant fraction of energy in 
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electrons and magnetic fields, relative to the thermal energy of the shocked particles, no 
choice of u m and u a (including the case when v m « u a ) can explain the rapid decay of the 
radio flux together with the observed optical afterglow light curve. Our calculations are thus 
in disagreement with the model proposed by Galama et al. (1999); essentially the same 
difficulty exists with the models used by Wang, Dai & Lu (1999) and Dai & Lu (1999). Shi 
& Gyuk (1999) overcome this difficulty by appealing to special geometry (rapid interaction 
of the blast wave with a cloudlet). Apart from invoking tremendous energy in the burst 
(10 55 erg) their model is unable to account for the smoothly decaying optical afterglow light 
curve. 

Galama et al. (1999) speculate that the broadening of the characteristic frequencies 
(v a , Vm-, v c ) may lead to a rapid decay of the radio flux. In our opinion, broadening the 
characteristic frequencies will smooth the predicted light curve but this does not help 
address the fundamental problem of rapid decay. In the framework of the forward shock 
model, the radio flare can be accounted only if the characteristic frequencies can be evolved 
much faster than that given by standard dynamics. In contrast, the reverse shock model 
provides a natural and consistent explanation for the radio flare. 



5. Radio Emission From the Forward Shock 

Accepting the conclusion that the radio flare arises in the reverse shock we now discuss 
the radio emission at late times, t 3> £*. At such late times, we expect some radio emission 
from the forward shock but with conventional properties, i.e., radio flux slowly rising to 
a maximum value (f m ) at time t m ~ weeks, and then decaying in much the same way as 
the optical afterglow. This is an inevitable consequence of the observed optical afterglow 
emission. 
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The highest measured afterglow flux (excluding the ROTSE observations that we 
associate with the reverse shock emission) is the Gunn r-band flux, f r = 100 /iJy at epoch 
t\ ~ 3.7 hr (Kulkarni et al. 1999). Clearly, f m > f r and t r , the epoch at which the r band 
flux peaks is < t\. Thus we expect a peak flux > f r in the 8.46-GHz band and this peaks 
at epoch t r (v r /v*) 2 ' 3 < 210 d; here, v T is the center frequency of the r-band, ~ 4 x 10 14 
Hz. From this it follows that we should have seen at least 44/xJy in < 40 d after the burst. 
As can be seen from Figure [I] and Table [TJ (also §|3|) our observations are sensitive enough 
to exclude such strong emission. Radio emission can be suppressed by invoking u a > v*. 
However, this condition would also suppress the radio emission from the reverse shock and 
thus be unable to account for the radio flare (see §||). 

Kulkarni et al. (1999) have interpreted the break in their r-band afterglow emission 
as due to a jet. The lack of significant radio emission is in excellent agreement with the 
jet hypothesis. Let tb be the time at which the jet's lateral dimensions become visible to 
the observer and/or the spreading of the jet begins. If v < v a then the radio flux is oc t° 
but oc t -1 / 3 , otherwise (Sari, Piran & Halpern 1999). We do note that these scalings for a 
spreading jet are different from those derived by Rhoads (1998) who assumed v m <C v a . 

From the r-band observations we obtain £& — 2.1 d (Kulkarni et al. 1999). Up 
to epoch tb, the radio flux increases as t 1 ' 2 and the expected radio flux at epoch tb is 
10(/ m /100AiJy) 13/9 /iJy- Limiting to t > t b we fit a t x / 3 power law to observed fluxes after 
subtracting the contribution from the reverse shock (see Figure [TJ also Sari & Piran 1999 
for details of the reverse shock model). The best fit model yields 22(t/4) -1 / 3 /iJy from 
which we infer f m = 170 /iJy. At the 95% confidence level, the radio emission is constrained 
to be less than 40(t/tfc)~ 1/ ' 3 /iJy which implies f m < 260 /iJy An additional and independent 
constraint comes from our late time (May 1999) observations (these observations were not 
used in the fits shown in Figure |l|). At the 95% level, our model predicts 11 /iJy which 
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should be compared to the observed mean flux of 7 ± 7 fiJy (see Table [1]). 

Combining the f m inferred from the radio observations with the r-band power law 
decay (flux oc t -1 ' 1 ), we find t r > 1.6 hr. Earlier in this section we noted some constraints 
on f r and t r . Taken together, we find 100 < f m < 260 /iJy and 1.6 < t r < 3.7 hr. Several 
authors (e.g. Galama et al. 1999, Castro-Tirado et al. 1999) have proposed that the 
ROTSE light curve smoothly blends into the optical afterglow emission. This hypothesis 
yields f m ~ 5 mJy, certainly inconsistent with our result above. In our opinion, the 
observational basis for this assertion is questionable given the differing slopes of the ROTSE 
light curve and the optical afterglow light curve. 

We end by noting that early time radio observations could potentially offer a unique 
insight into the physics of the reverse shock (Sari & Piran 1999, Meszaros & Rees 1999). 
Motivated thus, we have re-examined archival data on other bursts and found evidence for 
a radio flare in GRB 970828. The specific incidence of radio flares of 1:8 appears to be 
significantly larger than that of prompt optical emission obtained by ROTSE or LOTIS. 

DAF thanks R. Hjellming and Marc Verheijen for generously giving up their VLA time 
so that a fast response could be made. RS is supported by a Fairchild Fellowship. 
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Table 1. Radio Observations of GRB 990123 



Epoch (UT) 


Telescope 


t (min) 


Band (mm) 


S±a (//Jy) 


1999 Jan 23.63 


VLA 


15 


36 


62±32 


1999 Jan 24.65 


VLA 


15 


36 


260±32 


1999 Jan. 24.68 


JCMT 


47 


1.35 


-4640 ± 2410 


1999 Jan 26.58 


VLA 


10 


36 


53±39 


1999 Jan 27.62 


VLA 


26 


36 


18±25 


1999 Jan. 27.82 


JCMT 


37 


1.35 


0.00 ± 1700 


1999 Jan. 27.88 


JCMT 


75 


0.85 


-3100 ±1100 


1999 Jan 28.32 


VLA 


37 


36 


22±25 


1999 Jan. 29.88 


JCMT 


60 


0.85 


630 ± 1700 


1999 Jan 31.27 


VLA 


76 


36 


49±19 


1999 Feb 02.52 


VLA 


25 


36 


4±26 


1999 Feb 04.43 


VLA 


78 


36 


14±16 


1999 Feb. 04.81 


JCMT 


40 


0.85 


4790 ± 1440 


1999 Feb. 05.83 


JCMT 


90 


0.85 


700 ± 1200 


1999 Feb. 06.84 


JCMT 


45 


0.85 


1950 ± 1820 


1999 Feb 07.36 


VLA 


48 


36 


10±23 


1999 Feb. 07.83 


JCMT 


15 


0.85 


4100 ± 2900 


1999 Feb. 07.85 


JCMT 


30 


1.35 


-800 ± 2500 


1999 Feb 11.37 


VLA 


109 


36 


16±16 


1999 Feb 14.44 


VLA 


44 


36 


17±18 


1999 Feb 18.43 


VLA 


31 


36 


27±17 


1999 Feb 21.34 


VLA 


136 


36 


24±12 


1999 Feb 22.29 


VLA 


145 


36 


-6±11 
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Table 1 — Continued 



Epoch (UT) 


Telescope 


t (min) 


Band (mm) 




1999 Feb 25.37 


VLA 


46 


36 


-3±18 


1999 Feb 27.37 


VLA 


86 


36 


14±14 


1999 Mar 04.27 


VLA 


132 


36 


8±13 


1999 May 06.29 


VLA 


110 


36 


6±10 


1999 May 07.41 


VLA 


170 


36 


9±10 



Note. — The columns are (left to right), (1) Start of observations (UT date), (2) Telescope 
(3) Duration of observations, (4) Band of observations, (5) Mean flux density and rms (//Jy) 
at the position of the radio transient. 
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Fig. 1. — The observed and model lightcurves at 8.46 GHz. Triangles are upper limits being 
defined here as the measured flux +2cr. The weighted averages over ranges of days are shown 
by boxes; the vertical size is ±cr. (A) The heavy curves are the expected light curves from 
our model. The left curve is the predicted emission from the reverse shock model of Sari & 
Piran (1999); hereafter model "R". The right curve is the sum of the radio emission from 
the reverse and the forward shock advocated in the text (model "F", see also Figure 
but with a jet geometry. We assume that the jet begins to spread at epoch 2.1 d (see §||). 
(B) The light solid lines are the light curves computed for a conventional spherical forward 
shock model with parameters favored by Galama et al. (1999), v a ~ v m ~ 30 GHz (model 
"Fg", see §^j). The smooth curve is the curve obtained from exact calculations of this model 
whereas the other curve is built up from asymptotic formulae. (C) The curves with dashed 
lines refer to the F G model but with a jet geometry similar to that discussed in (A) above. 
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GRB 990123: epoch 1.2 day 
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Fig. 2. — The observed broad-band spectrum close to the epoch of Jan 24.65. The radio 
data are 4.84 GHz (Galama et. al. 1999, epoch of Jan 24.46), 8.46 GHz (this paper, Jan 
24.65), 15 GHz (Galama et al. 1999, Jan 25.26), 86 GHz and 232 GHz (Jan 25.18; Galama 
et al. 1999) and 222 GHz (this paper, Jan 24.68). Triangles are 3a upper limits. The optical 
point is obtained from the power law interpolation of the r-band light curve (index —1.1) 
of Kulkarni et al. (1999). The light solid curve labelled "Fq" is the expected forward shock 
contribution from model "Fq" (see caption to Figure 0). The spectral slope at high frequency 
(optical) is set to —0.75, consistent with the observed rate of decay of the optical afterglow 
(see Kulkarni et al. 1999). See the text for details on the computation of the curve. The 
heavy solid curve "R" is the predicted radio emission from reverse shock model (Sari & Piran 
1999). The heavy solid curve "F" is the radio afterglow emission from the forward shock; 
the parameters for this model are displayed in the Figure. 



